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Abstract 

This  document  presents  the  software  requirements  for  the  Network  Common  Data 
Format  (netC'DF)  interfaces  that  are  being  used  to  modularize  the  Experimental  Cen¬ 
ter  for  Mesoscale  Ocean  Prediction  (ECMOP)  being  built  at  the  Institute  for  Naval 
Oceanography  (INO).  ECMOP  is  currently  under  development.  It  has  been  designed 
to  be  a  bread-board  test  bed  for  ocean  model  evaluation  and  comparisons. 
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1  Introduction 


The  goal  of  this  document  is  to  define  the  netCDF  environment  within  which  the  Experi¬ 
mental  Center  for  Mesoscale  Ocean  Prediction  (ECMOP).  which  is  being  developed  at  the 
Institute  for  Naval  Oceanography  (INO).  will  operate.  Towards  that  end  we  present  a  data 
dictionary  defining  each  data  item  that  might,  appear  in  any  of  our  net  CDF  documents. 

We  will  consider  only  the  numerical  ocean  prediction  models  STEM  [1]  and  PR1DER 
[2][l].  and  the  data  assimilation  module  ASSM  [fi],  in  this  version  of  the  document.  It  is 
expected  that  the  document  will  be  appended  as  other  models  are  added  to  ECMOP.  There 
are  separate  sections  for  both  SPEM  and  PRIDER.  and  each  new  model  will  need  to  provide 
a  section  describing  its  own  net-CDF  variable  name  mapping,  also. 

We  suggest  another  use  for  this  document,  however.  New  code  for  t  he  models  can  adopt 
eit  her  this  exact  data  dictionary  or  a  logical  reduction  of  it.  In  such  a  situation  t  he  problems 
of  leading  a  new  code  would  be  eased  for  people  familiar  with  the  INO  data  dictionary  for 
ocean  data  items. 

The  netC'DF  interface  has  been  designed  to  allow  separate  distinct  task  environments 
to  communicate  through  a  “format  independent"  common  data  format  (CDF).  The  "net" 
indicates  that  these  files  can  be  read  on  any  computer  and  across  net  works. 

It  swms  appropriate  to  justify  the  need  for  this  document.  The  netCDF  lile  format 
is  format  independent,  and  carries  its  variable  names  with  it.  Reading  and  understanding 
the  file  contents  then  depends  on  the  variable  names  used  and  the  dimension  names  that 
define  t  he  shape  of  those  variables.  By  coupling  a  data  dictionary  with  a  netCDF  file,  the 
contents  of  that,  file  are  easier  to  understand  and  use  through  a  simple  interface,  wit  hout 
much,  assistance.  This  is  our  goal  here. 

2  Information  Description 

A  description  of  the  netCDF  can  be  found  in  [7].  INO's  purpose  in  moving  to  netCDF  is  to 
provide  a  more  llexible  environment  for  the  PM)2  assessment.  At  that  assessment  we  will  be 


t  est  ing  many  aspects  of  each  model.  After  the  assessment  we  need  to  he  able  to  combine  t  in* 
more  successful  aspects  of  each  model  into  a  composite  model  which,  hopefully,  is  greater 
than  the  sum  of  its  parts.  We  anticipate  that  the  netCDF  generality  will  make  any  redesign 
work  proceed  much  quicker. 

The  netCDF  software  comes  with  two  utilities,  negen  and  ncdunip.  These  can  be  used, 
in  a  limited  sense,  instead  of  a  universal  interpreter  for  the  CDF.  Given  a  file  descript  ion  in 
( 'l.)F,  ncfjt  n  will  generate  a  set  of  ('  or  FORTH  AN  subroutine  calls  that  write  I  lie  appropriate 
information  to  a  netCDF  lile.  The  INO  written  utility  inoconr  (Sec.  8.2)  will  change  the 
"write"  subroutine  calls  to  "read"  subroutine  calls.  This  utility  produces  a  corresponding 
set  of  subroutine  calls  that  read  the  netCDF  file. 

'There  is  a  one-to-one  mapping  between  the  netCDF  read  and  write  subroutine  calls,  and 
fortunately  the  parameters  are  the  same  for  each  corresponding  call.  The  utility  maps 

ncvptl  ->  ncvgtl 
ncvplc  ->  ncvglc 
nevpt  ->  nevgt 

But  other  subroutine  calls  must  be  ignored.  The  CDF  structure  provides  for  a  deline 
mode  or  state.  In  define  state  the  attributes  and  dimensions  are  established.  This  needs 
to  be  done  only  for  the  initial  write  routine.  Also  there  are  options  that  rename  attributes 
within  the  file  headers,  and  synchronize  (checkpoint  )  the  current  status  of  the  write  program 
with  the  file  contents.  We  do  not  allow  any  of  these  subroutines  to  be  called  in  the  read 
subroutine  interface  to  the  netCDF  file.  The  in  or  o  nr  utility  comments  out  these  calls.  The 
commented  out  routine  calls  are  listed  below. 


ncredf 

-  enter  define  mode 

ncendf 

-  exit  define  mode 

nesne 

-  force  a  checkpoint  to  the  file 

nccre 

-  create  a  netCDF  file 

neddff 

-  create  a  dimension 

ncapt 

-  create  an  attribute 

ncaptc 

-  create  am  attribute 

ncadel 

-  delete  am  attribute 

ncaren 

-  rename  am  attribute 

Two  subroutines  are  then  available  to  read  and  write  t  he  netC'DF  file  and  these  subrou¬ 
tines  rail  be  linked  and  then  loaded  along  with  the  routines  that  wish  to  access  this  data. 
Fa ch  model  then  needs  only  one  set  of  interface  routines,  instead  oi  a  separate  set  ol  routines 
for  each  file  format. 

Of  course  the  routine  must  be  regenerated  each  time  a  change  has  to  be  made  to  the 
tile  format,  but  the  utility  usage  is  straightforward  [7].  After  regenerating  the  netODF  calls 
using  the  provided  utilities,  they  must  be  relinked  into  the  routine  using  them. 

In  the  future  we  plan  to  add  the  capability  of  reading  and  writing  to  and  from  a  buffer 
instead  of  a  lile.  This  amounts  to  formatting  an  internal  buffer.  This  would  save  uu needed 
lile  access  time  when  running  on  very  large  supercomputers. 

We  now  present,  the  detailed  description  of  the  problem  the  netCDF  interfaces  are  be¬ 
ing  used  to  solve.  This  includes  human-machine  interaction,  and  hardware  and  software 
elements. 

r) 


Figure  2.1:  Full  Connect  ion  wit  hout  net  CDF 

2.1  Data  Flow  Diagrams 

YVc  are  using  the  common  data  format  to  reduce  an  0(\e)  complexity  problem  to  a  problem 
of  O(n).  With  separate  formats  for  each  (ile.  all  tasks  must  have  communication  jackets 
for  all  other  tasks.  For  example,  with  three  models  and  three  auxiliary  general  purpose 
processing  t  asks,  we  would  need  nine  jackets  to  convert  between  the  three  varieties  of  each 
component,  M  and  T  (Fig.  2.1).  If  we  were  to  add  two  more  varieties  of  models  and  tasks 
t  he  number  of  jackets  required  would  jump  to  25. 

You  can  see  that  this  is  a  fast  growing  communications  problem.  With  the  net  CDF  jackets 
the  situation  becomes  more  manageable.  With  a  universal  translator  the  communication 
paths  would  be  as  in  Fig.  (2.2). 

That  is,  full  connectivity  can  be  achieved  with  an  0{i>)  linear  number  of net  CDF  jackets. 
For  a  system  that  needs  full  connectivity  real-time  the  advantages  of  simple  connections 
far  out  weigh  the  cost  of  the  translations. 

For  the  current  INO  system.  ECMOP.  such  real-time  connectivity  is  not  required.  This 
is  fortunate  because  the  authors  could  not  find  an  existing  universal  translator,  and  had  no 
ambitions  to  write  their  own.  Instead  we  will  use  the  utilities  to  generate  only  the  specific 
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Figure  2.2:  Full  Connection  With  uetCDF 

nelODF  jacket  needed  for  the  system  we  are  configuring  for  a  test  .  For  each  run  we  will  use 
a  connection  as  in  Fig.  (2.3) 


Figure  2.3:  Load  Time  Connections  With  ECMOP 

This  nets  simplicity  of  connectivity  at  run-tmie  without  the  cost  involved  for  the  universal 
translator.  Yet  it  leaves  us  upwards  compatible  so  that,  if  and  when  a  universal  translator 
becomes  available  we  would  like  to  have  the  option  to  use  it. 

2.2  ECMOP  Overview 

The  Institute  for  Naval  Oceanography  (INO)  was  founded  in  1985.  by  a  Secretary  of  the  Navy 
Initiative,  to  accomplish  mesoscale  ocean  prediction  on  a  global  basis  and  to  incorporate 
university  advances  in  ocean  modeling  into  Naval  Oceanography.  This  document  is  part  of 
the  INO  Experimental  Center  for  Mesoscale  Ocean  Prediction  (ECMOP). 

ECMOP  consists  of  three  major  sub-systems:  the  Verification  Module  (YEHMOD).  the 
Data  Sub-System  (DASS).  and  the  Visualization  Module  (Y1SMOD).  The  data  assimilation 
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(ASSM)  routines  described  within  this  document  sit  between  the  DASS  (iml  tin-  mode!  code. 
Our  Ion!  model  code,  the  semi-spectral  primitive  equation  model  SPEM.  is  presented  in 
llcdstrom  [I], 

A  preprocessing  utility  routine  retrieves  the  necessary  data  from  the  DA.SS  through  EM¬ 
PRESS.  a  Data  Base  Management  System  (DBMS),  and  formats  it  for  ASSM.  ASSM  then 
ingests  the  data  real  time  and  transforms  if  to  the  requested  variables  and  format  needed 
to  feed  the  model  routine.  This  communication  takes  place  within  the  nct.CDF  environment 

[']• 

Fig,  ( 'J .  I )  shows  a  high  level  diagram  of  the  ECMOP  system.  The  outer  box.  which  is 
labeled  "iiel.CDF  Network  Environment",  serves  as  a  reminder  that  the  subsystem  column- 
nioations  will  take  place  within  the  netCDF  environment.  These  communication  links  are 
marked  in  the  diagram  with  arrows  denoting  the  directional  flow  of  the  paths  of  communica¬ 
tion.  Rectangular  boxes  are  used  to  indicate  the  major  software  components,  man-machine 
interactions,  storage  devices,  and  any  form  of  input  or  output. 

For  example.  Dr.  Louise  Perkins,  one  of  the  authors  of  this  document,  is  responsible  for 
Ibis  design  document,  the  software  characteristics  of  the  ASSM  module,  the  coding,  and  for 
the  unit  testing  within  the  ECMOP  system  at  1NO. 

Our  design  philosophy  was  to  design  and  define  in  detail  the  complete  functionality 
expected  from  the  individual  modules.  But  implementation  proceeded  iteratively,  beginning 
from  the  simpliest  functionality  and  evolving  towards  the  highest  levels. 

The  DBMS  is  EMPRESS  and  is  installed  on  a  SUN  SPARC  station  at  Stennis  Space 
Center  in  Mississippi.  The  Atmospheric  Sciences  Directorate  of  the  Naval  Oceanographic 
and  Atmospheric  Research  Laboratory  (NOAKL),  formerly  NEPRF.  has  been  developing  a 
DBMS  system  for  atmospheric  research  [3].  INO  systems  personnel  have  taken  maximum 
advantage  of  their  expertise  and  their  software,  and  are  using  their  DBMS  design.  A  prepro¬ 
cessor  extracts  ail  data  needed  for  data  assimilation  for  each  model  run  from  I  his  DBMS,  or 
from  a  data  file. 
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not('f)F  Network  Environment 


Current  timing  estimates  oil  the  Massachusetts  Institute  of  Technology  Supercomputer 
facility  (MITSF)  Cray-2  indicate  that  approximately  two  Cray  minutes  are  required  on  one 
processor  for  each  day  the  STEM  model  is  advanced.  The  Cray-2  at  the  Massachusetts 
Institute  of  Technology  (MIT)  is  not  currently  heavily  utilized  so  that  the  CPC  time  and 
the  wall  clock  time  are  comparable.  Hence  we  consider  this  a  real-t  ime  requirement  lor  our 
data  access. 

The  DBMS  resides  on  a  separate'  computer,  and  network  file  transport  is  needed  to  up¬ 
load  the  data  onto  the  larger  supercomputer.  For  this  reason  we  use  a  preprocessor  and 
assemble  the  net  CDF  data  liles  before  the  actual  run  begins.  The  indeliniteness  of  future 
data  assimilation  data  needs,  coupled  with  a  lack  of  timing  results  for  the  DBMS  and  the 
network  access  path,  made  this  our  salest  design  choice. 

Model  output  can  be  stored  by  the  DBMS,  using  a  postprocessor,  and  \  EH  MOD  will 
be  able  to  obtain  the  output  fields  and  compare  one  against  another.  Y  ISM  01.)  is  used  to 
display  graphical  results  after  or  during  a  test. 

3  High  Level  Design 

As  an  example,  we  describe  the  design  and  development  of  the  ASSM  module.  The  ASSM 
module  is  being  implemented  in  four  iterative  steps  (Fig.  (To))-  The  lirst  functional  assim¬ 
ilation  program,  the  innermost  square,  is  called  ASSM  HC  (Hard  Cotied).  It  provides  the 
fund  ioualily  to  incorporate  current-  meter  data  available  at  numerical  grid  points.  ASSM  IK  ' 
was  a  proof-of-concept .  cotied  by  Dr.  Roberta  Young.  It  was  hardcoded  without  flexibility, 
and  no  documentation  is  available.  It  is  not,  considered  part  of  ECMOP. 

Coding  for  the  ASSM  BS  (Basic  STEM)  was  the  second  iteration  and  has  been  used 
exclusively  with  the  STEM  model. 

Next  the  connection  to  a  preprocessor  that  retrieves  data  was  implemented.  This  is  the 
third  level  of  iteration.  ASSM  DB. 
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Next  the  netCDF  interfaces  for  AS8M  were  implemented.  This  completed  the  ASSM 
utility:  the  fourth  iteration  is  referred  to  as  ASSM  CDF. 

A  unit  test  of  STEM  within  the  ECMOP  environment  will  he  performed  when  t  he  FA 
MOP  system  is  ready.  This  is  expected  to  locate  the  interface  bugs. 

ASSM  CDF  can  be  expanded  to  support,  multiple  models.  Additional  users  are  expected 
to  be  the  Princeton  Model  (PEDAM),  configured  for  the  North  Atlantic,  and  the  Data. 
Assimilation  Research  and  Transition  (DART)  program  of  NOARL.  These  groups  will  enter 
FCMOP  with  verification  coding  and  interface  standards  completed:  ASSM  will  provident! 
interface  for  nudging  data  assimilation  into  their  codes. 

Future  iterations  may  be  funded  for  transition  to  the  Naval  Oceanography  Command. 
These  would  be  general  purpose  versions  which  r  ;|uire  more  rigorous  testing  and  validation 
than  is  described  here. 

4  Data  Dictionary 

'Hits  section  presents  the  data  dictionary  we  are  using.  All  words  are  abbreviated  to  rum¬ 
or  loss  let  ters.  Less  than  four  is  acceptable  only  for  words  less  than  four  let  ters,  or  for  well 
established  commonly  accepted  abbreviations.  The  abbreviations  selected  followed  some 
general  guidelines.  The  first  guideline  was  to  drop  all  vowels  when  possible,  except  when  the 
first  letter  is  a  vowel.  We  rarely  drop  the  first  letter  in  the  word  unless  it  is  a  number  or  a 
commonly  used  English  language  abbreviation.  If  we  still  need  more  reduction,  we  next  pick 
the  most  important  consonants  to  retain.  Third  a  visual  inspection  for  readability  was  used. 

The  shorter  abbreviations  are  three  or  less  letters.  It  is  possible  to  have  a  hierarchy  of 
shorter  abbreviations.  W'c  first  tried  to  shorten  an  abbreviation  by  removing  letters  from  the 
end  of  the  long  abbreviation.  When  this  was  not  unique,  we  retained  the  dominant  sounds 
of  the  accented  (stressed)  syllables  [0]. 


The  type  field  is  keyed: 
desc  -  description 
attr  -  attribute 
rvar  -  real  variable 
ivar  -  integer  variable 
cvar  -  character  variable 
bvar  -  bit  or  logical  variable 
gvar  -  generic  variable 
op  -  operation 
unit  -  unit  of  measure 


Full.Name 

Abbr_Name 

Short 

Type 

abbreviated 

abbr 

abr , ab 

desc 

absolute 

abs 

al 

op 

access 

acc 

as 

op 

acoustic 

acst 

ac 

attr 

action 

actn 

atn.an 

op 

active 

actv 

av 

desc 

advection 

advc 

ad,  a 

rvar 

ad  j  o  int 

adjt 

aj 

op 

age 

age 

attr 

agency 

agcy 

age 

desc 

air 

air 

ai 

rvar 

alt imetry 

altm 

al 

attr 

amount 

amt 

am 

desc 

anemometer 

anmt 

anm 

attr 

l’J 


angle 

ang 

ag 

op 

anomaly 

anora 

desc 

ansi 

ansi 

desc 

area 

area 

ar 

op 

array 

arry 

ara 

desc 

ascending 

asc 

attr 

assimilation 

assm 

ass ,  as 

op 

associative 

assc 

act 

op 

atmosphere 

atms 

af 

attr 

attached 

att 

desc 

attribute 

attr 

at 

attr 

avhrr 

avhr 

attr  Adv .  Very  Hi  Resol .  Raidometer 

background 

bkgr 

attr 

band 

band 

attr 

bandwidth 

bdwh 

desc 

baroclinic 

bcnc 

be 

attr 

barotropic 

btpc 

bt 

attr 

base 

base 

desc 

beaufort 

bfrt 

attr  Beaufort  wind  strength 

becomes 

bcom 

op 

begin 

bgn 

b 

desc 

biology 

biol 

desc 

bit 

bit 

desc 

bitstream 

bstr 

bvar 

bogus 

bgus 

bgs 

desc 

boresight 

bors 

attr  sensor  orientation 

bottom 

bot 

bm 

desc 

boundary 

bdry 

bd 

desc  edge  or  limits 

It 


brightness 

brgt 

rvar  luminance 

byte 

byte 

desc 

calculate 

calc 

op 

calls 

call 

op 

callsign 

csgn 

attr  call-sign  of  station  or  ship 

calm 

calm 

desc 

catalog 

cat 

desc  descriptive  list  of  options 

category 

catg 

desc  division  of  classification 

cell 

cell 

desc 

Celsius 

cels 

ce 

attr 

center 

cntr 

desc 

change 

chg 

desc  modification  or  alteration 

channel 

chan 

attr 

characteristic 

chr 

desc 

charge 

chrg 

cirrus 

cirr 

ci 

desc  type  of  high  cloud 

classification 

cl  as 

desc  systematic  organization 

clear 

clr 

desc 

climatology 

clmo 

cmo 

attr  data  pertaining  to  climate 

cloud 

cld 

desc 

cm 

cm 

attr  central  memory 

coast 

cst 

attr 

code 

code 

cd 

op  an  encoding 

coefficient 

coef 

c 

gvar 

column 

col 

desc 

command 

cmd 

op  order  or  instruction 

comment 

cmt 

desc 

common 


com 


desc 


compilation 

cmpl 

op 

complexity 

cplx 

op 

component 

cmp 

attr 

computer 

cptr 

desc 

concentration 

cone 

attr 

condition 

cond 

desc  restricting  or  modifying  factor 

conditional 

cndl 

op 

confidence 

enf  i 

attr  degree  of  certainty 

constant 

enst 

attr 

contains 

entn 

op 

cont inuous 

cont 

desc 

contour 

entr 

desc  isoline 

control 

Ctrl 

desc 

convective 

conv 

rvec 

convergence 

envg 

op 

conversion 

envn 

op  transformation  between  states 

copies 

copy 

op 

Coriolis 

corl 

cor, cl  op 

correction 

corr 

op  to  meet  a  standard 

count 

ent 

attr  Numeric  value  for  cardinality 

country 

ctry 

desc 

coverage 

cvrg 

desc  extent  of  coverage 

cp 

cp 

desc  central  processor 

creates 

ere 

op 

creation 

ertn 

op 

criticality 

crit 

desc 

crossection 

xset 

op 

crossreference 

xref 

op 
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cube 

cube 

cb 

desc 

current 

curr 

rvar 

curviture 

curv 

rvar 

cycle 

eye 

ivar 

dampen 

damp 

op 

data 

data 

d 

gvar 

date 

date 

dt 

attr  time  in  days 

day 

day 

attr 

decay 

decy 

dey 

deck 

deck 

desc  of  punched  cards 

decode 

ded 

op 

deep 

deep 

rvar 

delayed 

dlyd 

op 

density 

dens 

dn 

rvar 

depth 

dpth 

dph.dp 

rvar 

depression 

dprs 

rvar  a  curviture  upwards 

descending 

dsd 

op 

describes 

dscr 

desc 

description 

dsc 

desc  Textual  appendage 

destination 

dest 

desc 

device 

dev 

desc 

dewpoint 

dwpt 

rvar 

diameter 

dia 

rvar 

dictionary 

diet 

desc 

difference 

dif 

df 

op 

diffusivity 

diff 

di 

op 

dimension 

dim 

dm 

desc 

direction 

dir 

desc 
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directory 

drct 

desc 

disk 

dsk 

desc 

displacement 

dspl 

rvar 

disposition 

disp 

desc 

distance 

dis 

ds 

rvar 

distribution 

dist 

desc 

diurnal 

drnl 

desc 

domain 

dom 

desc  limits  of  variable 

Dopier 

dplr 

desc 

east 

e 

desc 

echo 

echo 

desc 

eddy 

eddy 

ey 

desc 

effect ive 

eff 

desc 

elapsed 

elps 

desc 

elevation 

elev 

elv , ev 

rvar 

end 

end 

ed 

desc 

energy 

engy 

rvar 

executable 

xqt 

xq 

desc 

executes 

exec 

op 

expanded 

xpnd 

desc 

external 

xtnl 

xt 

desc 

factor 

fctr 

ftr.fr 

attr 

fine 

fine 

desc 

first 

1 

desc 

flag 

flag 

fl 

desc  boolean  flag 

float 

f  lot 

fit 

desc 

flux 

flux 

f 

rvar 

formal 

frml 

desc 

IS 


friction 

fret 

fre 

rvar 

front 

frnt 

desc 

full 

full 

desc 

function 

fuac 

fun.fc 

op 

geographic 

geog 

desc  taxonomy  of  surface  of  earth 

geosat 

gsat 

desc 

get 

get 

g* 

op 

GHz 

ghz 

gh 

attr 

global 

glbl 

desc  worldwide,  comprehensive 

gravity- 

grvt 

gvt 

rvar 

grid 

grid 

gr  >g 

desc 

ground 

grnd 

desc  surface  of  earth 

g+' 

g* 

op  greater  than 

guess 

gues 

op  estimation 

heavy 

hvy 

desc 

heat 

heat 

he 

rvar 

height 

ht 

rvar 

hertz 

hz 

h 

attr 

high 

hi 

desc 

horizontal 

horz 

desc 

hour 

hour 

rvar 

hydrography 

hydr 

attr 

identity 

id 

desc  ident.  for  object  or  attribute 

image 

imge 

im 

desc 

index 

indx 

ix 

desc 

inactive 

inac 

desc 

influence 

inf  1 

if 

desc 

information 


info 


desc 


infrared 

inf  r 

attr 

ingest 

ngst 

op 

input 

iput 

ip 

op 

insertion 

isrt 

in 

op 

instaneous 

inst 

desc 

instruction 

inst 

itr 

op 

intern  ad 

itnl 

it 

desc 

intervad 

itvl 

ivl.il, 

i  desc 

inverted 

inv 

iv ,  v 

op 

ir 

ir 

attr 

is 

is 

desc 

isentropic 

isnt 

desc 

isobaric 

isbr 

desc 

kilogram 

kilo 

kg 

unit 

kinematic 

knmt 

desc 

Lagrangiaui 

lagn 

attr 

land 

land 

lnd 

desc 

last 

last 

desc 

latent 

ltnt 

desc 

latitude 

lat 

rvar 

left 

left 

1 

desc 

length 

lnth 

gvar 

level 

lvl 

lv 

desc 

light 

lit 

desc 

list 

list 

desc 

lit 

lit 

attr 

load 

load 

desc 

location 

lctn 

desc 

infra-red 

exists 

equal  potential  temperature 


any  list  of  words  or  numerals 
info  about  lit  data 
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desc 


logical 

logc 

lg 

long 

long 

longitude 

Ion 

low 

low 

lw 

It 

It 

lu 

lu 

maps 

maps 

mps 

marine 

marn 

mask 

mask 

msk.ms 

mass 

mass 

matches 

mtch 

maximum 

max 

mean 

mean 

mn 

meander 

mand 

method 

mthd 

mt 

meters 

mtrs 

mtr.mr 

middle 

mid 

m 

millibar 

mb 

minimum 

min 

minutes 

mnts 

missing 

miss 

mixed 

mixd 

model 

modi 

mdl 

modifies 

mdfy 

molecular 

mlcr 

momentum 

mmtm 

mm 

moored 

mord 

nadir 

ndr 

rvar 

desc 

op  less  than 
desc  logical  unit 
op 

desc 

op 

rvar 

op 

op 

op 

desc 

attr 

unit 

desc 

attr 

op 

attr 

desc 

desc 

desc 

op 

attr 

rvar 

attr 

desc  point  below  surface  observation 
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name 

name 

desc 

literal  value  of  object 

next 

next 

desc 

north 

n 

desc 

0  degrees 

northeast 

ne 

desc 

45  degrees 

northwest 

nw 

desc 

315  degrees 

not 

not 

op 

nudging 

nudg 

ndg.ng 

op 

number 

num 

nu 

ivar 

index  into  series  of  objects 

numeric 

nmrc 

nine 

desc 

observation 

ob 

0 

desc 

occupies 

occp 

desc 

ocean 

ocn 

desc 

oe 

oe 

attr 

info  about  orbital  element 

operational 

oper 

desc 

operator 

op 

attr 

operation  to  perform 

ordinal 

ord 

desc 

output 

oput 

op 

op 

particle 

prtc 

pc 

desc 

past 

past 

desc 

path 

path 

pth,ph 

desc 

per 

per 

pe 

op 

period 

perd 

desc 

duration  of  interval  event 

permanent 

perm 

desc 

pk 

pk 

desc 

acronym  -  primary  key 

point 

pt 

desc 

points 

pts 

P 

desc 

precipitable 

prcl 

desc 

present 

prst 

desc 

pressure 

pres 

rvar 

prevailing 

prvl 

Pr 

desc 

previous 

prev 

pv 

desc 

primary 

prim 

desc 

print 

prnt 

prt 

op 

probability 

prbl 

op 

profiles 

prof 

prf ,pf 

desc 

put 

put 

op 

q 

q 

desc  terrain  resolution 

qualitative 

qltv 

desc 

radiation 

rdtn 

rvar 

radius 

rad 

ra 

rvar 

Raf  os 

raf  s 

raf 

desc 

ramp 

ramp 

op 

random 

radm 

desc 

range 

rang 

desc 

rate 

rate 

rat 

desc 

read 

read 

rd 

op 

references 

ref 

desc 

relative 

rltv 

desc 

remark 

rmrk 

desc  addendum  for  clarity 

report 

rpt 

desc 

request 

rqst 

rq 

op 

requester 

reqr 

desc 

requires 

req 

desc 

research 

rsch 

desc 

resembles 

rsmb 

desc 

resolution 

rslt 

desc 

responsible 

resp 

desc 

return 

rtn 

rt 

op 

right 

rht 

r 

desc 

Rossby 

rsby 

desc 

salinity 

salt 

sal , si 

rvar 

satellite 

sat 

st 

attr 

save 

save 

sv 

op 

scale 

scl 

sc 

op 

sea 

sea 

se 

desc 

second 

2 

desc 

secondary 

scdy 

desc 

seconds 

sec 

desc 

sensible 

snsb 

desc 

sensor 

sens 

desc 

shallow 

shlw 

desc 

shear 

shr 

rvar 

ship 

ship 

desc 

sigma 

sgma 

rvar 

signif icant 

sig 

desc 

simple 

smpl 

sp 

desc 

size 

size 

sz 

attr 

smooth 

smth 

desc 

sounding 

sndg 

desc 

solar 

sol 

desc 

sonic 

snc 

desc 

sounder 

sndr 

desc 

south 

s 

desc  180  degrees 

southeast 

se 

desc  135  degrees 

southwest 

sw 

desc  225  degrees 

space 

spc 

desc 

spawns 

spwn 

op 

speed 

spd 

rvar 

squared 

sqrd 

sq 

op 

start 

strt 

desc 

station 

stn 

desc 

step 

step 

desc 

stream 

strm 

desc 

stress 

strs 

str 

rvar 

subset 

sbst 

ss 

desc 

subjective 

sbj 

desc 

subordinate 

sub 

desc 

sum 

sum 

op 

superior 

supr 

desc 

surface 

surf 

sf 

desc 

Sverdrup 

sver 

attr 

temperature 

temp 

tp 

rvar 

temporary 

tmp 

desc 

terrain 

terr 

desc 

text 

text 

txt 

desc  qualitative  info,  in  text 

thermocline 

thcl 

desc 

third 

3 

desc 

threshold 

thld 

desc 

time 

tm 

t 

desc  absolute  time 

topography 

topo 

desc 

total 

totl 

ttl 

desc 

track 

trek 

tk 

desc 
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transformation 

trfm 

op 

transient 

trns 

desc 

transport 

tprt 

desc 

turbulence 

trbl 

trb ,tl 

desc 

type 

type 

ty 

attr  selection  from  a  category 

u 

vldm 

rvar  velocity  in  first  dimension 

ucomponent 

ucmp 

desc 

unexpended 

unxp 

desc 

unit 

unit 

u 

desc 

units 

unts 

US 

desc 

update 

updt 

ud 

op 

upper 

uppr 

up 

desc 

urban 

urbn 

urb 

desc 

uses 

use 

desc 

usn 

usn 

desc  US  Navy 

V 

v2dm 

rvar  velocity  in  second  dimension 

valid 

vald 

desc 

value 

val 

desc  numeric  represent,  of  quantity 

variable 

vrbl 

vb 

desc 

variance 

var 

desc 

vcomponent 

vcmp 

desc 

velocity 

vel 

v 

rvar 

vertical 

vert 

vr 

desc 

viscosity 

vise 

vs 

rvar 

visible 

vis 

desc 

vorticity 

vort 

rvar 

w 

v3dm 

rvar  velocity  in  third  dimension 

warming 

wrmg 

desc 
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water 


wtr 


desc 


wave 

wav 

desc 

wban 

wban 

desc  weather  bureau  Air  Force  Navy 

weather 

wthr 

desc 

weight 

wght 

wt 

desc 

west 

w 

desc  270  degrees 

wet bulb 

wetb 

attr  wetbulb  thermometer 

whitecap 

wcap 

desc  crest  of  wave  is  foam 

width 

wdth 

desc 

wind 

wind 

wnd 

rvar 

wmo 

wmo 

desc  World  Meteoro.  Organization 

word 

word 

desc 

write 

wrt 

wr 

op 

xtom 

xtom 

desc  computation  excludes  toms 

yaw 

yaw 

op  rotation  about  vertical  axis 

year 

year 

yr 

desc 

zenith 

zen 

desc  vertically  above  observation 

zone 

zone 

desc 

zonal 

znl 

desc 

5  net  CDF  Jackets 

The  ucK.’DF  jackets  are  responsible  for  mapping  model  variable  names  to  t  he  net  CDF  long 
and  short  abbreviated  names  given  in  the  data  dictionary  in  Sec.  (I).  Complex  names, 
consisting  of  underscore-concatenated  long  form  abbreviations,  are  used  for  the  net  CDF  lile 
variable  names,  whenever  possible.  Short  form  abbreviations  are  concatenated  without  the 
underscore.  In  lengt  h  limited  programming  languages  the  shorter  abbreviations  are  used  in 
the  netCDF  lile. 


27 


This  section  includes  the  currently  defined  mappings  for  both  the  SPEM  and  Princeton 
models  as  well  as  for  the  ASSM  module.  All  can  be  viewed  as  components  in  the  Et.’MOP 
environment. 

5.1  SPEM  to  netCDF  Mapping 

We  map  the  SPEM  variable  names,  which  are  listed  in  the  column  “model",  to  the  “Short- 
form"  column  abbreviations  available  in  the  neK'DF  file  because  SPEM  is  a  EOUTllAN 
program. 


long 

Model 

netCDF 

dt 

time_step 

dtyp 

data_type 

disp 

flag.dcd 

h 

ht 

ipre 

flag_prev 

isave 

f lag_3ave 

iupd 

flag_updt 

idat 

date 

irw 

flag.rd.wrt 

ire 

f lag_rtn_code 

img 

data_iflc 

irpt 

logc.unit.rpt 

ixb 

arry.size_bgn_l.dim 

ixe 

arry.size_end_l.dim 

ivb 

arry.3ize_bgn_2.dim 

iye 

arry.3ize_end_2.dim 

short 

netCDF  Description 


dt 

dty 

fldcd 

ht 

flpr 

f  lsv 

flud 

date 

flrdwr 

f lrtcd 

diflc 

lurpt 

bldm 

edldm 

b2dm 

ed2dm 
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izb 

arry_size_bgn_3_dim 

b3dm 

ize 

arry _ s ize_end_3_dim 

ed3dm 

iuptr 

l_dim_vel_ob_lctn_l_dim 

vllodm 

ivptr 

2_dim_vel_ob_lctn_l_dim 

v21odm 

irptr 

dens_ob_lctn_l_dim 

dnlodm 

juptr 

l_dim_vel_ob_lctn_2_dim 

vl2odm 

jvptr 

2_dim_vel_ob_lctn_2_dim 

v22odm 

jrptr 

dens_ob_lctn_2_dim 

dn2odm 

kuptr 

l_dim_vel_ob_lctn_3_dim 

vl3odm 

kvptr 

2_dim_vel_ob_lctn_3_dim 

v23odm 

krptr 

dens_ob_lctn_3_dim 

dn3odm 

loca 

lctn 

lctn 

1 

l_dim_grid_size 

ldmgr 

m 

2_dim_grid_size 

2dmgr 

mtppt 

max_topo_size 

toposz 

numu 

l_dim_vel_ob_size 

vlodjn 

numv 

2_dim_vel_ob_size 

v2odm 

numr 

dens  _ob_  size 

dnodm 

n 

3_dim_grid_size 

gr3dm 

nptlg 

num_pts_lon 

nplon 

nptlt 

num_pts_lat 

nplat 

psi 

vort 

vort 

rampuv 

per 

per 

resol 

rslt 

rslt 

rho 

dens 

dens 

rlllt 

low_left_lat 

lvllat 

rllng 

low_left_lon 

lwllon 

robs 

dens.ob 

dno 

grid_size 

grid_size 

period  to  interpolate  thru 


coverage,  geographic_location 
coverage,  geographic.location 
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ros2 

cnst_rsby_rad 

rsbyrd 

Rossby  radius 

rclm 

dens.clmo 

dnclmo 

rdcy 

arry_data_if lc 

aradif 

rdata 

dens_data 

dnd 

rudec 

1  _d  im_  vel  _c  o  ef  _  as  sm 

vlcas 

rvdec 

2_dim_vel_coef _assm 

v2cas 

rrdec 

dens_coef _assm 

dncas 

scalf 

scl_fctr 

sclftr 

ssbathy  sbst_topo 

sstopo 

sslllt 

sbst_low_left_lat 

sslwl2 

domain_values  or 

coverage 

ssllg 

sbst_low_left_lon 

sslwll 

domain_values  or 

coverage 

ssurlt 

sb s t _ uppr_ rht _ 1 at 

ssupr2 

domain.values  or 

coverage 

ssurlg 

sbst_uppr_rht_lon 

ssuprl 

domain_values  or 

coverage 

t 

curr_tm 

currtm 

tobs 

ob_tm 

obtm 

tdata 

elps.tm 

elpstm 

time  since  last 

data  assimilation 

u 

l_dim_vel 

vldm 

uclm 

l_dim_vel_clmo 

vldmcl 

udata 

l_dim_vel_data 

vldmd 

urlat 

uppr_rht_lat 

uprlat 

urlng 

uppr_rht_lon 

uprlon 

uobs 

l_dim_vel_ob 

vldob 

vobs 

2_dim_vel_ob 

v2dob 

V 

2_dim_vel 

v2dm 

vdata 

2_dim_vel_data 

v2dmd 

vclm 

2_dira_vel_clmo 

v2dmcl 

w 

3_dim_vel 

v3dm 
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5.2  PRIDER  to  netCDF  Mapping 

The  mapping  for  the  PRIDER  package  is  given  in  this  subsection.  We  map  the  PRIDER 
variable  names,  which  are  listed  in  the  column  “Model",  to  the  Short  form  abbreviation 
because  PRIDER  is  a  FORTRAN  routine. 


Model 

netCDF 

Short  Description 

aam2d 

eddy_vi.sc 

eddyvs 

alat 

grid_cntr_lat 

grlat 

alon 

grid_cntr_lon 

grlon 

ang 

ang_l_dim_2_dim 

angl2 

Angle  between  1  and  2  dimensions 

art 

area_l_dim_2_dim_grid 

areal2 

Area  defined  by  cells  in  1  and  2  dim 

aru 

area_cntr_vel_ l_dim 

vldm 

arv 

area_cntr_vel_2_dim 

v2dm 

a  am 

knmt.visc 

knmtvs 

Kinematic  Viscosity 

cbc 

coef_bot_frct 

cbtfrc 

cor4 

corl 

corl 

curv42d 

sum_corl 

sumcrl 

covrhn 

grvt _ rdt n_n _bdry 

gvtnbd 

Gravity  radiation  at  northern  boundary 

covrhe 

grvt  _  rdtn_e_bdry 

gvtebd 

d 

ht 

ht 

dtef 

vert_prof 

vertpf 

dum 

l_dim_vel_mask 

vldmms 

dvm 

2_dim_vel_mask 

v2dmms 

dx 

l_dim_dis_invl 

dsldm 

dy 

2_dim_dis_invl 

ds2dm 

dt 

3_dim_dis_time_step 

ds3dmt 

dte 

xtnl_time_step_sec 

xtnlt 

dti 

itnl_time_step_sec 

itnlt 

dz 

3_dim_dis_invl 

ds3dmi 

dzz 

3„dim_dis_mid 

ds3dmm 

dzr 

3_dim_dis_inv_invl 

ds3drav 

egf 

elev 

elev 

egb 

elev_prev_t ime_  st ep 

elevts 

el 

elev_xtnl 

elevxt 

elb 

elev_prev_tiae_step_xtnl 

elevpv 

et 

elev_int 

elvint 

etb 

elev_prev_t ime_step_itnl 

elvitt 

etf 

elev_next_time_step_itnl 

elvit 

eln 

sea_lvl_n_bdry 

selvnb 

ele 

ssa_lvl_e_bdry 

selveb 

els 

sea_lvl_s_bdry 

selvsb 

fluxua 

advc_l_diro 

aldm 

fluxva 

advc_2_dim 

a2dm 

f  sm 

sea_lvl_mask 

selvms 

ftc 

damp_temp_salt_f ctr 

tpsal 

grav 

grvt 

grvt 

h 

3_dim_depth 

dph3dm 

i 

l_dim_grid_indx 

gldmix 

im 

l_dim_grid_size 

gldmsz 

int 

ixtl_time_step 

ixtlt 

iprint 

prnt_invl 

prtiv 

jm 

2_dim_grid_size 

g2dmsz 

j 

2_dlm_grid_indx 

g2dmix 

k 

3_dim_grid_indx 

g3dmix 

kh 

vert.diff _coef 

vrdic 

km 

vert  Jaunt 

vrknmt 

kq 

vert_trbl 

vrtrbl 

1 

trbl_lnth_scl 

trbl 

mode 

type.calc 

calcty 

nts 

next_t ime.step 

nextt 

pts 

prev_time_step 

prevt 

psi 

strm_func 

strraf c 

q2 

2_trbl_engy 

trbl2 

q2b 

2_trbl_engy_pts 

trb2pt 

q21 

2_trbl_engy_lnth_scl 

trb2sc 

q21b 

2_trbl_engy_lnth_scl_pts 

tlscpt 

rho 

dens 

dens 

rhof 

dens_pts 

denspt 

rmean 

area_mean 

areamn 

ramp 

ramp_strt_fctr 

rajnpfr 

3 

salt 

salt 

sb 

salt_pts 

saltpt 

sbe 

salt.ejadry 

salebd 

sbn 

salt_n_bdry 

salnbd 

sbs 

salt_s_bdry 

salsbd 

sbw 

salt_w_bdry 

salwbd 

smean 

mean_salt 

mnsalt 

t 

temp_cels 

tempce 

tb 

temp_cels_pts 

tpcept 

tmean 

temp_cels_mean 

tpcemn 

tprnu 

wght.fctr 

wtf ctr 

time 

time 

time 

tbe 

temp_e_bdry 

tpebd 

tbn 

temp_n_bdry 

tpnbd 

tbs 

temp_s_bdry 

tpsbd 

tbw 

temp_w_bdry 

tpwbd 

uabe 

btrp_vel_e_bdry 

btvebd 

u 

l_dim_vel 

vldm 

uf 

l_dim_vel_pts 

vldmpt 

ub 

l_dim_vel_num_time_step 

vldmnu 

ua 

l_dim_vel 

vldm 

uab 

l_dim_vel_pts 

vldmpt 

uaf 

l_dim_vel_num_time_step 

vltm 

umol 

mlcr_diff 

mlcrdi 

va 

2_dim_vel 

va2dm 

V 

2_dim_vel 

v2dm 

vf 

2_dim_vel_pts 

v2dmpt 

vb 

2_dim_vel_num_time_step 

v2dmnu 

vh 

vert.prof 

vertpf 

vhp 

vert.prof _pts 

vrpfpt 

vabn 

btrp_vel_n_bdry 

btvnbd 

vabs 

btrp_vel_s_bdry 

btvsbd 

vab 

2_dim_vel_pts 

v2dmpt 

vaf 

2  _  d im_ vel _num_ t ime_  s t ep 

v2dmnu 

wssurf 

salt_flux_surf 

salfsf 

wtsurf 

heat_flux_surf 

htfsf 

wubot 

l.dim_mmtm_f lux.bot 

mmlbt 

wvbot 

2_dim_mmtm_flux_bot 

mm2bt 

wusurf 

l_dim_ramtm_flux_surf 

mmlsf 

Momentum  flux  at  bottom  in  1  dimension 
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wvsurf  2_dim_mmtm_flux_surf  mm2sf 

w  3_dim_vel  v3dm 

z  3_dim_dis  ds3dm 

zz  3_dim_dis_mid_invl  ds3iv 

5.3  ASSM  Routine 

The  mapping  for  the  ASSM  package  is  given  in  this  subsection.  We  map  the  ASSM 


names,  which  are  listed  in  the  column 
ASSM  is  a  FORTRAN  routine. 


Model 

netCDF 

Long  Name 

asfun 

assm_func 

biold 

biol.data 

elm 

clmo_arra 

data 

data 

dens 

dens 

dt 

time.step.size 

gridpt 

grid.pt 

idat 

curr.date 

inf ogt 

inf o.get 

inst 

inst 

io 

iput.oput.flag 

ipre 

prev.flag 

ire 

rtn.code.flag 

irng 

rang. if lc 

“Model'*,  to  the  Sliortl'orm  abbreviation 


netCDF 

Shortform  Description 


asfun  what  to  do 

biold 

cmara 

data 

dens 

tstpsz 

gridpt 

currdt 

inf ogt 

inst 

ipopfl 

prevfl 

rtncd 

rangif 


variable 

because 


is 

intr 

isave 

save.flag 

iupd 

updt_flag 

1 

size_l_dim 

lagnd 

lagn_flot_info 

lctn 

lctn 

m 

size_2_dim 

mdlasd 

modl_vrbl_assm_data 

mdldsc 

modl_grid_dsc 

mdlvb 

modl_vrbl_val 

mean 

mean 

model 

modi 

n 

size_3_dim 

nudclm 

nudg_clmo 

nump 

num_pts 

ob 

ob 

pres 

pres 

rdcy 

rate_decy 

rho 

dens_arra 

spce 

spce 

sver 

tprt.sver 

t 

time 

temp 

temp 

time 

time 

tobs 

time.ob 

u 

vel_l_dira_arra 

V 

vel_2_dim_arra 

vel 

vel 

intr 

savef 1 

updtfl 

szldm 

lagnd 

lctn 

sz2dm 

mdlasd 

mdldsc  what  it  is 

mdlvr 

mean 

mdl  who  to  assimilate  for 

sz3dm 

nudgcm 

nupts 

ob 

pres 

ratdcy 

dens 

spce 

sver 

t 

temp 

time 

tob 

vldm 

v2dm 

vel 
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V 


vel_3_dim_arra 


v3dm 


wndstr  wnd.strs  wndstr 

5.3.1  Data  Space 

There  are  many  “standard"  data  structures  in  data  space.  The  best  structure  for  each  data 
type  with  respect  to  storage  needs,  retrieval  and  update  needs  are  still  open  issues.  In 
this  section,  with  an  eve  towards  the  future,  we  set  down  the  contents  of  groups  of  data 
by  collection  instruments  or  commonly  associated  data  types.  The  instrument  is  listed  and 
then  some  of  the  data  types  we  might  expect  to  be  associated  with  that  instrument  are  listed 
under  it . 


ADP: 

acst_dplr_prof  (adp) 
vel.prof _vert  (vprfvr) 

Gross  Feature  Locations  (Bogus  Maps) : 
mndr 
eddy 

fmt.trak  (frnttk) 

axis 

n_wall 


Hydrography : 

Hydr 

sea_surf_ht  (sesfht) 
sea_surf _terap  (sesftp) 
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Inverted  Echo  Sounders : 
inv_echo_sndr  (ies) 
thcl.dpth 

Path  Characteristics: 
path_char  (pthchr) 
dspl 

ang_prtc_path  (agpcph) 
curvature  (curv) 

Satellite  Altimetry: 
sat.altm  (satal) 

sea_surf_ht  (sesfht) 
infr  (infr) 


6  Constant  Values 

For  actual  numbers,  the  value  can  be  concatenated  with  the  unit  or  type.  For  ('sample 
1(1  —  miniitf- s.  or  (lie  associated  short  form  10 mnt9.  can  be  used  lor  a  constant  time  value 
of  ten  minutes. 

7  Constraints  on  Data  Items 

Some  words  are  not  allowed  because  other  words  are  equivalent  and  "better  in  some  sense. 
The  following  words  are  not  to  be  added  to  the  data  dictionary  in  Sec.  (2): 

altitude  -  use  height  or  elevation 


average  use  mean 


datetime  use  time 


definition  use  description 
increment  use  interval 

mode  use  description  or  type 

quantity  use  amount 

8  Functional  Description 

This  section  describes  the  procedural  details  of  using  nelCDF  within  EC.' MOP.  There  it  re 
several  options.  Binding  the  netC'DF  interface  to  the  model  can  he  made  at  compile  time 
by  generating  new  subroutines  for  any  changes.  Binding  can  occur  at  load  time  by  select  ing 
one  of  several  subroutines  to  use,  and  binding  can  occur  at  run  time  by  making  dynamic 
decisions  about  the  contents  of  a  nelCDF  lile.  For  the  FY92  time  frame,  we  suggest  load 
time  binding  using  makefiles  within  the  UNIX  environment.  This  is  the  most  eost  effect  ive 
in  t  he  short  run. 

8.1  Processing  Narrative 

This  section  will  walk  us  through  the  algorithms  available  for  use  within  each  utility  and 
each  subroutine. 

The  nelCDF  provides  several  data  formats: 

•  byte 

•  character 

•  short 

•  long  ( 3*2  bit) 

•  lloal  (32  bit  IEEE  floating-point) 

•  double  (64  bit) 
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The  INO  ECMOP  environment  lias  defniod  some  coimnou. attributes: 


•  units 

•  long- name 

•  valid- range 

•  valid-mill 

•  valid- max 

•  sca  le-  factor 


•  add-offset 

•  missing-value 

•  ('-format 

-  F OUTRAN-  format 

•  title 

•  history 

Given  the  contents  of  data  files  that  exist  or  that  need  to  be  created,  the  user  must 
genera  te  a  header  fde  that  describes  the  contents  of  the  file.  The  Sl’EM  netCDF  hist  or  .  file 
interface  is  written  with  the  following  FORTRAN  program.  The  corresponding  ODL  file  is 
presented  after  the  program. 
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************************************************** 


*  Name:  MCSSTgencdf  * 

*  Programmer:  Phuong  Vaui  Leach  * 

*  Company:  Institute  for  Naval  Oceanography  * 

*  Date:  July  19,  1991  * 

************************************************** 

*  * 

*  This  FORTRAN  program  creates  a  netCDF  file  * 

*  in  which  the  scientific  data  cam  be  accessed  * 

*  and  shared  in  a  form  that  is  self-describing  * 

*  and  network-transparent .  * 

*  * 

************************************************** 

program  hist 


*  links  to  netCDF  libratry 

include  'netcdf.inc' 

********************************** 

*  declares  variables  for  parameter 

integer  nhist 
integer  charlen 
integer  num.long 
integer  num.lat 
integer  nura_ levels 
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*  defines  maximum  dimensions 

parameter  (nhist  =  4) 
parameter  (charlen  =  36) 
parameter  (num_long  =  130) 
parameter  (num.lat  =  130) 
parameter  (num_levels  *5) 

*  declares  hypercube 

integer  strt2D(2) ,  cnt2D(2) 
integer  strt3D(3) ,  cnt3D(3) 

*  declares  hypercube  data 

data  strt2D  /l,  1/ 

data  strt3D  /l,  1,  1/ 

data  cnt2D  /num.long,  num_lat/ 

data  cnt3D  /num_long,  num.lat,  num_ levels/ 

*  declares  return  error  code 

integer  iret 

*  declares  netCDF  id 

integer  cdfid 

*  delcares  dimension  ids 

integer  nhistdim,  charlendim 
integer  lpdira,  mpdim,  npdim,  mdim 
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*  declares  variable  ids 


integer  xiuid,xivid,xitid,xisid,psiid,identid 

integer  xlid , elid ,hid ,f id ,uvnu2id 

integer  snu2id , tnu2id ,uvnu4id , snu4id ,tnu4id , rhoOid 

integer  gid , rdrgid .hminid ,hmaxid , rf acid , rdmpid , pmid 

integer  pnid , dndxid , dmdeid , dhdxid , dhdeid , dxid , deid 

integer  xpid , ypid ,piid, sigid , cpid , cf id , cdid , cdzid 

integer  cint id , csumid , iprmid , f prmid , mgopt id 

integer  xminid,yminid,xmaxid,ymaxid 

integer  sclmid,tclmid,uclmid,vclmid 

integer  dt id, time id 


*  declares  variable  shapes 
integer  diml(l) 
integer  dim2(2) 
integer  dim3(3) 


♦**%** ***************%***♦*****♦******* ********** 


*  declares  data  variables 

real  xiu(num_long ,num_lat ,nura_levels) 
real  xiv(num_long,num_lat ,nun_levels) 
real  xit(num_long,num_lat ,num_levels) 
real  xis(num_long,num_lat,num_levels) 
real  psi(num_long,num_lat) 
character  ident (charlen) 
real  xl 
real  el 
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real  h(num_long,num_lat) 

real  f (num_long,num_lat) 

real  uvnu2 

real  snu2 

real  tnu2 

real  uvnu4 

real  snu4 

real  tnu4 

real  rhoO 

real  g 

real  rdrg 

real  hmin 

real  hmax 

real  rfac 

real  rdmp 

real  pm 

real  pn 

real  dndx 

real  dmde 

real  dhdx 

real  dhde 

real  dx 

real  de 

real  xp 

real  yp 

real  pi 

real  sig 

real  cp 
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real  cf 
real  cd 
real  cdz 
real  cint 
real  csum 
integer  iprm 
integer  fprm 
integer  mgopt 
real  xmin 
real  ymin 
real  xmax 
real  ymax 
real  sclm 
real  tclm 
real  uclm 
real  vclm 
real  dt 
real  time 

4^************>M************************************************ 

*  This  routine  generates  a  netcdf  file  from  a  data  file.  * 

**************************************************************** 

*  creates  netcdf  file/  enters  define  mode 

cdfid  *  nccre( '/epoch/data3/leach/hist .cdf ' ,  ncclob,  iret) 


*  defines  dimensions 
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nhistim  *  ncddef (cdf id,  ’nhist',  nhist,  iret) 
charlendim  *  ncddef (cdf id,  'charlen’,  charlen,  iret) 
lpdim  ■  ncddef (cdfid,  ’num.long’,  num_long,  iret) 
mpdim  *  ncddef (cdfid,  ’num_lat',  num.lat,  iret) 
npdim  *  ncddef (cdfid,  'num.levels' ,  num.levels,  iret) 

*  defines  dbms  variables 
dim3(3)  3  npdim 
dim3(2)  *  mpdim 
dim3(l)  3  lpdim 

xiuid  *  ncvdef (cdfid, ’xiu’ .NCFLOAT, 3, dim3, iret) 
xivid  3  ncvdef (cdfid, ’xiv’ .NCFLOAT, 3, dim3, iret) 
xitid  =  ncvdef (cdf id, ’xit ’ .NCFLOAT, 3, dim3, iret) 
xisid  *  ncvdef (cdf id, ’xis’ .NCFLOAT, 3, dim3, iret) 

dim2(2)  3  mpdim 
dim2(l)  3  lpdim 

psiid  3  ncvdef (cdfid, 'psi’ .NCFLOAT, 2, dim2, iret) 
diml(l)  3  charlendim 

identid  3  ncvdef (cdfid, ' ident ' ,NCCHAR, l,diml , iret) 

xlid  3  ncvdef (cdfid, ’xl’ , NCFLOAT, 0,0, iret) 
elid  3  ncvdef (cdfid, ' el’ .NCFLOAT, 0,0, iret) 


dim2(2)  3  mpdim 
dim2(l)  3  lpdim 
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hid  *  ncvdef (cdf id, 'h' , NCFLOAT, 2, dim2, iret) 
fid  *  ncvdef (cdfid, 'f' .NCFLOAT, 2, dim2, iret) 

uvnu2id  *  ncvdef (cdfid, ’uvnu2' , NCFLOAT, 0,0, iret) 
snu2id  =  ncvdef (cdfid, ’snu2’ .NCFLOAT, 0,0, iret) 
tnu2id  =  ncvdef (cdf id, ’tnu2’ .NCFLOAT, 0,0, iret) 
uvnu4id  *  ncvdef (cdfid, 'uvnu4' .NCFLOAT, 0,0, iret) 
snu4id  =  ncvdef (cdfid, ’snu4' .NCFLOAT, 0,0, iret) 
tnu4id  *  ncvdef (cdfid, ’tnu4’ .NCFLOAT, 0,0, iret) 
rhoOid  =  ncvdef (cdf id, 'rhoO ’, NCFLOAT, 0 ,0 , iret) 
gid  *  ncvdef (cdf id, 'g' .NCFLOAT, 0,0, iret) 
rdrgid  *  ncvdef (cdfid, ’rdrg' .NCFLOAT, 0,0, iret) 
hminid  *  ncvdef (cdfid, ’hmin’ .NCFLOAT, 0,0, iret) 
hmax id  »  ncvdef ( cdf id , ’ hmax ’ .NCFLOAT ,0,0, iret ) 
rfacid  ■  ncvdef (cdf id, 'rfac' .NCFLOAT, 0,0, iret) 
rdmpid  *  ncvdef (cdfid, 'rdmp' .NCFLOAT, 0,0, iret) 
pmid  *  ncvdef (cdfid, ’pm' .NCFLOAT, 0,0, iret) 
pnid  =»  ncvdef  (cdfid, 'pn' .NCFLOAT,  0,0,  iret) 
dndxid  ■  ncvdef (cdfid, 'dndx' .NCFLOAT, 0,0, iret) 
dmdeid  =  ncvdef (cdf id, 'dmde' , NCFLOAT, 0 ,0 , iret) 
dhdxid  =  ncvdef (cdfid, 'dhdx' .NCFLOAT, 0,0, iret) 
dhdeid  =  ncvdef (cdfid, 'dhde' .NCFLOAT, 0,0, iret) 
dxid  *  ncvdef ( cdf id, 'dx’, NCFLOAT, 0,0, iret) 
deid  *  ncvdef (cdfid, ’de’ .NCFLOAT, 0,0, iret) 
xpid  ■  ncvdef (cdfid, ’xp’ .NCFLOAT, 0,0, iret) 
ypid  *  ncvdef (cdfid, 'yp' .NCFLOAT, 0,0, iret) 
piid  ■  ncvdef (cdfid, 'pi' .NCFLOAT, 0,0, iret) 
sigid  *  ncvdef (cdfid, 'sig' .NCFLOAT, 0,0, iret) 


cpid  *  ncvdef (cdf id, 'cp' .NCFLOAT, 0,0, iret) 
cfid  *  ncvdef (cdfid, ’cf’ .NCFLOAT, 0,0, iret) 
cdid  *  ncvdef (cdfid, 'cd' .NCFLOAT, 0,0, iret) 
cdzid  ■  ncvdef (cdf id, ' cdz’ .NCFLOAT, 0,0, iret) 
cintid  ■  ncvdef (cdf id, ’cint’ ,NCFLOAT, 0,0, iret) 
csumid  *  ncvdef (cdfid, ’csum’ .NCFLOAT, 0,0, iret) 
iprmid  =  ncvdef (cdf id, ’iprm’ .MCLONG, 0,0, iret) 
fprmid  »  ncvdef (cdf id, ’fprm’ ,NCL0NG, 0,0, iret) 
mgoptid  *  ncvdef (cdfid, ’mgopt’ ,NCL0NG, 0,0, iret) 
xminid  =  ncvdef (cdf id, ’xmin’ .NCFLOAT, 0,0, iret) 
yminid  =  ncvdef (cdf id, 'ymin' ,NCFL0AT,0 ,0 , iret) 
xmaxid  =  ncvdef (cdfid, ’xmax' .NCFLOAT, 0,0, iret) 
ymax id  *  ncvdef ( cdf id , ’ ymax ' , NCFLOAT , 0 , 0 , ir et ) 
sclmid  =  ncvdef (cdfid, 'sclra' .NCFLOAT, 0,0, iret) 
tclmid  =  ncvdef (cdfid, ’tclm’ .NCFLOAT, 0,0, iret) 
uclmid  =  ncvdef  (cdfid,  ’uclm’  .NCFLOAT,  0,0,  iret) 
vclmid  *  ncvdef (cdf id, ’vclm’ .NCFLOAT, 0,0, iret) 
dtid  ■  ncvdef (cdf id, ’dt’ .NCFLOAT, 0,0, iret) 
timeid  *  ncvdef (cdf id, ’time’ .NCFLOAT, 0,0, iret) 

*  assigns  attributes  to  dbms  variables 

call  ncaptc(cdfid,xiuid, 'long.name' ,NCCHAR,31 , 

1  ’velocity  in  longitude  dimension’,  iret) 
call  ncaptc(cdf id.xivid, ’long.name’ ,NCCHAR,30, 

1  ’velocity  in  latitude  dimension’ ,  iret) 
call  ncaptc(cdf id.xitid, ’long.name’ .NCCHAR, 11 , 

1  'temperature' ,  iret) 

call  ncaptc(cdf id.xisid, ’long.name’ .NCCHAR, 11, 
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1  'temperature' ,  iret) 

call  ncaptc(cdf id,psiid, ’long_natme' ,NCCHAR,36, 

1  'horizontal  transport  stream  function',  iret) 
call  ncaptc(cdfid,identid, 'long.name' ,NCCHAR,40, 

1  'run  identification  (YYMMDD.type  of  data)',  iret) 
call  ncaptc(cdf id, slid, 'long_name' ,NCCHAR,25, 

1  'longitude  interval  length',  iret) 
call  ncaptc(cdf id, el id, 'long_name' ,NCCHAR,24, 

1  'latitude  interval  length',  iret) 
call  ncaptc(cdfid,hid, ’long_name’ ,NCCHAR,18, 

1  'bottom  depth  model',  iret) 
call  ncaptc(cdf id,f id, 'long_name' ,NCCHAR,15, 

1  'coriollis  force',  iret) 

call  ncaptc(cdf id,uvnu2id, 'long_name' ,NCCHAR,42, 

1  'laPlacian  coefficient  of  friction  velocity',  iret) 
call  ncaptc(cdfid,snu2id, 'long.name' ,NCCHAR,42, 

1  'laPlacian  coefficient  of  friction  salinity',  iret) 
cadi  ncaptc(cdfid,tnu2id, 'long.name' ,NCCHAR,45, 

1  'laPlacian  coefficient  of  friction  temperature', 

1  iret) 

call  ncaptc(cdf id,uvnu4id, 'long.name' ,NCCHAR,43, 

1  'biharmonic  coefficient  of  friction  velocity' , 

1  iret) 

call  ncaptc(cdfid,snu4id, 'long_name' ,NCCHAR,43, 

1  'biharmonic  coefficient  of  friction  salinity' , 

1  iret) 

call  ncaptc(cdf id,tnu4id, 'long_name' ,NCCHAR,46, 

1  'bihanmonic  coefficient  of  friction  temperature', 
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1  iret) 

call  ncaptc(cdf id.rhoOid, ’long_name’ ,NCCHAR,16, 

1  ’constant  density’,  iret) 
call  ncaptc(cdf id,gid , ’long_name’ ,NCCHAR,7, 

1  'gravity’,  iret) 

call  ncaptc(cdf id.rdrgid, ’long_name’ ,NCCHAR,23, 

1  ’bottom  drag  coefficient’,  iret) 
call  ncaptc(cdf id,hminid, ’long_name’ ,NCCHAR,20, 

1  'minimum  bottom  depth’,  iret) 
call  ncaptc(cdfid,hmaxid, ’long.name’ ,NCCHAR,20, 

1  'maximum  bottom  depth’,  iret) 
call  ncaptc(cdfid,rfacid, ’long_name' ,NCCHAR,38, 

1  ’simulation  conversion  factor  (.5  g/p.)',  iret) 
call  ncaptc(cdf id.rdmpid, ’long.name’ ,NCCHAR,27, 

1  'weights  for  nudging  forcing’ ,  iret) 
call  ncaptc(cdf id,pmid, 'longjaame’ ,NCCHAR,48, 

1  ’Jacobian  of  grid  transf ormation  in  Ion  dimension’, 

1  iret) 

call  ncaptc(cdf id.pnid, ’long.name’ ,NCCHAR,48, 

1  'Jacobian  of  grid  transf ormation  in  lat  dimension’, 

1  iret) 

call  ncaptc(cdfid,dndxid, ’long_name’ ,NCCHAR,43, 

1  ’Hermetian  of  "pn"  with  respect  to  longitude’,  iret) 
call  ncaptc(cdf id.dmdeid, ’long_name’ ,NCCHAR,42, 

1  'Hermetian  of  "pm"  with  respect  to  latitude’ ,  iret) 
call  ncaptc(cdf id.dhdxid, ’long_name’ ,NCCHAR,50, 

1  ’derivating  of  topography  with  respect  to  longitude’, 
1  iret) 


call  ncaptc(cdf id,dhdeid, 'long_name' ,NCCHAR,49, 

1  'derivating  of  topography  with  respect  to  latitude 
1  iret) 

call  ncaptc(cdf id,dxid, ’ long_name' ,NCCHAR,31 , 

1  'interval  in  longitude  dimension',  iret) 
call  ncaptc(cdfid,deid, 'long_name' ,NCCHAR,30, 

1  'interval  in  latitude  dimension',  iret) 
call  ncaptc(cdf id.xpid, 'long.name' ,NCCHAR,21, 

1  'longitude  grid  points',  iret) 
call  ncaptc(cdf id.ypid, 'long.name' ,NCCHAR,20, 

1  'latitude  grid  points',  iret) 
call  ncaptc(cdf id.piid, ' long_name' ,NCCHAR,2 , 

1  'pi',  iret) 

call  ncaptc(cdfid,sigid, 'long.name' ,NCCHAR,39, 

1  'stretched  coordinate  collocation  points’,  iret) 
call  ncaptc (cdf id ,cpid, 'long.name' ,NCCHAR,20, 

1  'polynomial  basis  set’,  iret) 
call  ncaptc(cdf id,cf id, ’long_name' ,NCCHAR,48, 

1  'matrix  of  basis  functions  at  the  sigma  locations' 
1  iret) 

call  ncaptc(cdf id,cdid, 'long_name' ,NCCHAR,44, 

1  'derivative  of  basis  with  respect  to  vertical’, 

1  iret) 

call  ncaptc(cdf id.cdzid, ’long_name' ,NCCHAR,44, 

1  'differentiation  of  a  model  field  in  vertical', 

1  iret) 

call  ncaptc(cdf id .cintid, 'long.name' ,NCCHAR,40, 

1  'integration  of  a  model  field  in  vertical',  iret) 


call  ncaptc(cdf id.csumid, ’long.name’ ,NCCHAR,24, 

1  ’Collocation  point  values’,  iret) 
call  ncaptc(cdf id.iprmid, ’long.name’ ,NCCHAR,23, 

1  ’mud2  integer  parameters ’ ,  iret) 
call  ncaptc(cdf id.fprmid, ’long_name' .NCCHAR, 21, 

1  ’mud2  float  parameters’ ,  iret) 
call  ncaptc(cdf id.mgoptid, ’long_name’ ,NCCHAR,32, 

1  'initialization  options  in  "init"’ ,  iret) 
call  ncaptc(cdf id.xminid, ’long_name' .NCCHAR, 17 , 

1  ’minimum  longitude’ ,  iret) 
call  ncaptc(cdf id ,yminid , ’ long.name ’ .NCCHAR, 16, 

1  'minimum  latitude',  iret) 
call  ncaptc(cdf id.xmaxid, ’long_name’ .NCCHAR, 17 , 

1  'maximum  longitude’,  iret) 
call  ncaptc(cdf id.ymaxid, ’long_name’ .NCCHAR, 16, 

1  'maximum  latitude’,  iret) 
call  ncaptc(cdf id.sclmid, ’long_name' .NCCHAR, 20, 

1  'salinity  climatology’ ,  iret) 
call  ncaptc(cdf id , tclmid , ’long_name’ .NCCHAR, 23, 

1  ’temperature  climatology’ ,  iret) 
call  ncaptc(cdf id.uclmid, ’long.name’ .NCCHAR, 43, 

1  ’velocity  in  longitude  direction  climatology’,  iret 
call  ncaptc(cdf id.vclmid, ’long_name’ , NCCHAR, 42, 

1  'velocity  in  latitude  direction  climatology',  iret) 
call  ncaptc(cdf id.dtid, ’long_name' .NCCHAR, 13, 

1  ’time  interval’,  iret) 
call  ncaptc(cdf id.timeid, ’long_name’ .NCCHAR, 4, 

1  ’time',  iret) 


*  leaves  define  mode 

call  ncendf (cdfid,  iret) 

*  initializes  variables 

do  30  i  *  l,lp 
do  20  j  *  l,mp 
do  10  k  *  l,np 
xiu(i, j ,k)  =0.0 
xiv(i,j ,k)  =0.0 
xit(i, j ,k)  =  0.0 
xis(i, j ,k)  =0.0 
10  continue 

psi(i.j)  =  0.0 
h(i  ,  j)  =  0.0 
f(i,j)  =  0.0 
20  cont inue 

30  continue 

do  40  i  =  l.charlen 
ident(i)  =  '  ’ 

40  continue 

11  =  0 
ml  =  0 
nl  =  0 
xl  =  0.0 


el  =  0.0 


uvnu2  *  0.0 
snu2  =0.0 
tnu2  =0.0 
uvnu4  =  0.0 
snu4  =0.0 
tnu4  =0.0 
rhoO  =0.0 
g  ■  0.0 
rdrg  =0.0 
hmin  =0.0 
hmax  =  0.0 
rfac  =0.0 
rdmp  =0.0 
pm  =  0.0 
pn  =  0.0 
dndx  =0.0 
dmde  =0.0 
dhdx  =0.0 
dhde  =0.0 
dx  =  0.0 
de  =  0.0 
xp  =  0.0 
yp  =  0.0 
pi  =  0.0 
sig  =  0.0 
cp  =  0.0 
cf  =  0.0 
cd  =  0.0 


cdz  *  0.0 


cint  *0.0 
csum  *  0.0 
iprm  =  0 
fprra  =  0 
mgopt  *  0 
xmin  =0.0 
ymin  =0.0 
xmax  =0.0 
ymax  =0.0 
sclm  =  0.0 
tclm  =  0.0 
uclm  =0.0 
vclm  =  0.0 
dt  =  0.0 
time  =0.0 

*  writes  values  into  netCDF  variables 

call  ncvpt (cdf id, xiuid,strt3D,cnt3D,xiu, iret) 

call  ncvpt(cdf id,xivid,strt3D,cnt3D,xiv, iret) 

call  ncvpt (cdf id , xit id , strt3D , cnt3D ,  xit , iret ) 

call  ncvpt (cdf id , xisid , strt3D , cnt3D ,  xis , iret ) 

call  ncvpt (cdf id, psi id, strt2D, cnt2D.ps i , iret) 

cadi  ncvptc (cdf id , ident id , 1 , charlen , ident , charlen , iret ) 

call  ncvptl(cdf id, Hid, 0,11, iret) 

call  ncvpt 1 (cdf id, ml id, 0, ml, iret) 

cadi  ncvptl(cdfid,nlid,0,nl,iret) 

call  ncvptl(cdfid,xlid,0,xl,iret) 


call  ncvpt 1 (cdf id, elid,0, el, iret) 
call  ncvpt (cdf id , hid , st rt 2D , cnt2D , h , iret ) 
call  ncvpt (cdf id,f id , strt2D , cnt2D ,f , iret) 
caill  ncvpt  1  (cdf id ,  uvnu2id ,  0 ,  uvnu2 ,  iret ) 
call  ncvpt 1 (cdf id , snu2id , 0 , snu2 , iret ) 
call  ncvpt 1 (cdf id , tnu2id , 0 , tnu2 , iret ) 
call  ncvpt 1 (cdf id ,uvnu4id , 0 , uvnu4 , iret ) 
call  ncvpt 1 (cdf id , snu4id , 0 , snu4 , iret ) 
call  ncvpt 1 (cdf id , tnu4id , 0 , tnu4 , iret ) 
call  ncvpt 1 (cdf id , rhoOid , 0 , rhoO , iret ) 
call  ncvpt 1 (cdf id, gid,0,g, iret) 
call  ncvpt 1 (cdf id , rdrgid , 0 , rdrg , iret ) 
call  ncvpt 1 (cdf id , hminid , 0 , hmin , iret ) 
call  ncvpt 1 (cdf id, hmaxid,0,hmax, iret) 
call  ncvpt 1 (cdf id, rf acid, O.rfac, iret) 
call  ncvpt 1 (cdf id , rdmpid , 0 , rdmp , iret ) 
call  ncvpt l(cdf id, pmid,0, pm, iret) 
call  ncvptl(cdfid,pn,0,pn,iret) 
call  ncvpt 1 (cdf id , dndxid , 0 , dndx , iret ) 
call  ncvptl(cdf id,dmdeid,0,dmde, iret) 
call  ncvptl(cdfid,dxid,0,dx,iret) 
call  ncvptl(cdf id,deid,0 ,deid,iret) 
call  ncvpt 1 (cdf id, xpid,0,xp, iret) 
call  ncvpt 1 (cdf id, ypid,0,yp, iret) 
call  ncvptl(cdfid,piid,0,pi,iret) 
call  ncvpt 1 (cdf id,sigid,0,sig, iret ) 
call  ncvpt 1 (cdf id , cpid , 0 , cp , iret ) 
cadi  ncvptl(cdfid,cfid,0,cf ,iret) 


call  ncvptl(cdfid,cdid,0,cd,iret) 
call  ncvpt 1 (cdf id , cdzid , 0 , cdz , iret) 
call  ncvpt 1 (cdf id , cintid , 0 , cint , iret ) 
cadi  ncvpt 1 (cdf id , csumid , 0 ,  csum , iret ) 
call  ncvpt 1 (cdf id , iprmid , 0 , iprm , iret ) 
call  ncvpt 1 (cdf id , f prmid , 0 , f prm , iret ) 
call  ncvptl(cdfid,mgoptid,0,mgopt ,iret) 
call  ncvpt 1 (cdf id , xminid , 0 , xmin , iret ) 
call  ncvpt  1  (cdf id , xmaucid ,  0 ,  xmaux ,  iret ) 
call  ncvpt  1  (cdf  id ,  ymaxid ,  0 ,  ymaLX ,  iret ) 
call  ncvpt 1 (cdf id , sclmid , 0 , sclm , iret) 
call  ncvpt 1 (cdf id , tclmid , 0 , t elm , iret ) 
call  ncvpt 1 (cdf id , uclmid , 0 , uclm , iret ) 
call  ncvpt 1 (cdf id , vclmid , 0 , vclm , iret ) 
call  ncvptl(cdfid,dtid,0,dt,iret) 
call  ncvpt 1 (cdf id , t imeid , 0 , t ime , iret ) 

*  closes  netCDF  file 

call  ncclos(cdf id,  iret) 

*  end  of  program 

end 

The  corresponding  CDL  is: 

net cdf  hist  { 
dimensions: 
nhist  =*  4  ; 


charlen  *  36  ; 
num_long  «  130  ; 
num_lat  *  130  ; 
num.levels  *  5  ; 

variables : 

float  xiu(num_levels,  num_lat,  num.long)  ; 

xiu : long_name  a  "velocity  in  longitude  dimension"  ; 

float  xiv(num_levels,  num_lat,  num_long)  ; 

xiv : long.name  ■  "velocity  in  latitude  dimension"  ; 

float  xit(num_levels,  num_lat,  num_long)  ; 

xit :  long_name  *  "temperature"  ; 

float  xis(num_levels ,  num.lat,  num_long)  ; 

xis : long_name  *  "temperature"  ; 

float  psi(num_lat,  num_long)  ; 

psi : long_name  *  "horizontal  transport  stream  function"  ; 
char  ident (charlen)  ; 

ident : long_name  *  "run  identification  (YYMMDD.type  of  data)"  ; 
float  xl  ; 

xl:long_name  *  "longitude  interval  length"  ; 
float  el  ; 

el:long_name  =  "latitude  interval  length"  ; 
float  h(num_lat,  num.long)  ; 
h:long_name  =  "bottom  depth  model"  ; 
float  f(num_lat,  num_long)  ; 
f:long_name  *  "coriollis  force"  ; 
float  uvnu2  ; 

uvnu2 : long_narae  *  "laPlacian  coefficient  of  friction  velocity" 
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float  snu2  ; 

snu2 : long.name  *  "laPlacian  coefficient  of  friction  salinity"  ; 
float  tnu2  ; 

tnu2:long_name  =  "laPlacian  coefficient  of  friction  temperature" 
float  uvnu4  ; 

uvnu4 : long_name  =  "biharmonic  coefficient  of  friction  velocity" 
float  snu4  ; 

snu4 : long_name  =  "biharmonic  coefficient  of  friction  salinity"  ; 
float  tnu4  ; 

tnu4 : long_name  =  "biharmonic  coefficient  of  friction  temperature 
float  rhoO  ; 

rhoO : long_name  =  "constant  density"  ; 
float  g  ; 

g:long_name  *  "gravity"  ; 
float  rdrg  ; 

rdrg:long_name  =  "bottom  drag  coefficient"  ; 
float  hmin  ; 

hmin:long_name  *  "minimum  bottom  depth"  ; 
float  hmax  ; 

hmax : long_name  *  "maximum  bottom  depth"  ; 
float  rfac  ; 

rfac:long_name  *  "simulation  conversion  factor  (.5  g/p.)"  ; 
float  rdmp  ; 

rdmp : long_name  =  "weights  for  nudging  forcing"  ; 
float  pm  ; 

pm:long_name  *  "Jacobian  of  grid  transf ormation  in  Ion  dimension' 
float  pn  ; 

pn:long_name  =  "Jacobian  of  grid  transformation  in  lat  dimension 


float  dndx  ; 

dndx:long_name  =  "Hermetian  of  \"pn\"  with  respect  to  longitude"  ; 
float  dmde  ; 

dmde : long_name  *  "Hermetian  of  \"pm\"  with  respect  to  latitude"  ; 
float  dhdx  ; 

dhdx : long.name  *  "derivating  of  topography  with  respect  to  longitude" 
float  dhde  ; 

dhde : long_name  =  "derivating  of  topography  with  respect  to  latitude" 
float  dx  ; 

dx: long_name  =  "interval  in  longitude  dimension"  ; 
float  de  ; 

de: long^name  =  "interval  in  latitude  dimension"  ; 
float  xp  ; 

xp:long_name  =  "longitude  grid  points"  ; 
float  yp  ; 

yp:long_name  *  "latitude  grid  points"  ; 
float  pi  ; 

pi:long_name  =  "pi"  ; 
float  sig  ; 

sig: long.name  *  "stretched  coordinate  collocation  points"  ; 
float  cp  ; 

cp : long_name  *  "polynomial  basis  set"  ; 
float  cf  ; 

cf:long_name  =  "matrix  of  basis  functions  at  the  sigma  locations"  ; 
float  cd  ; 

cd:long_name  =  "derivative  of  basis  with  respect  to  vertical"  ; 
float  cdz  ; 

cdz : long_name  *  "differentiation  of  a  model  field  in  vertical"  ; 
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float  cint  ; 

cint :long_name  *  "integration  of  a  model  field  in  vertical"  ; 
float  csum  ; 

csum : long.name  =  "Collocation  point  values"  ; 
long  iprm  ; 

iprm : long.name  *  "mud2  integer  parameters"  ; 
long  fprm  ; 

f prm : long.name  =  "mud2  float  parameters"  ; 
long  mgopt  ; 

mgopt : long_name  =  "initialization  options  in  V'initV’"  ; 
float  xmin  ; 

xmin:long_name  =  "minimum  longitude"  ; 
float  ymin  ; 

ymin : long_name  *  "minimum  latitude"  ; 
float  xraax  ; 

xmax : long.name  =  "maximum  longitude"  ; 
float  ymax  ; 

ymax:long_name  =  "maximum  latitude"  ; 
float  sclm  ; 

sclm:long_name  *  "salinity  climatology"  ; 
float  tclm  ; 

tclm:long_name  ■  "temperature  climatology"  ; 
float  uclm  ; 

uclm : long.name  »  "velocity  in  longitude  direction  climatology" 
float  vclm  ; 

vclm:long_name  *  "velocity  in  latitude  direction  climatology" 
float  dt  ; 

dt:long_narae  *  "time  interval"  ; 


float  time  ; 

time:long_name  =  "time"  ; 

} 


To  move  back  anti  forth  between  the  CDL  and  the  program,  the  user  moves  to  1  lie 
location  of  the  netCDF  utilities.  There  they  execute  the  negen  utility  with  options  to  create 
a  FORTRAN  program  that  reads  the  file  and  a  netCDF  file  itself  (if  desired).  The  utility 
call 

negen  -f  -o  <f ilename. write> 

generates  the  corresponding  FORTRAN  program.  A  C  language  option  is  also  available. 

8.2  INOCONV 

The  utility  INOCONV  (INO  Conversion)  can  now  be  run  on  this  FORTRAN  program  that 
writes  the  netCDF  file  to  generate  a  FORTRAN  program  that  reads  the  same  netCDF  file. 
This  is  executed  using  the  command  line: 

inoconv  <  routine. write. f  >  routine. read. f 

The  utility  program  is  straightforward  and  is  presented  here  in  its  entirety. 

c 


program  pttogt 
c  fortran  header 

c  module  name:  pttogt  (change  netedf  writes  (PuTs)  TO  reads  (GeTs)) 

c  module  type:  program  system  utility 

cmodule  environment:  kept  in  eemop  directory  for  general  use 

cmodule  description:  maps  nevp  to  nevg  calls  in  netedf  and  comments  out 

c  the  define  state 

eparameters  -  input :  none 


c  output:  none 

c  called  by:  none 

c  call:  none 

c  commons  used:  none 

cfiles  -  input:  *.in 

c  output:  *.out 

c  fortran  header 

character*l  linein(80) , lineout(80) 
logical  mcont 

open ( 18, file®' text . in’ , status®' old’ ) 
open (17, file®' text .out ’ .status® 'unknown’ ) 
mcont®. false, 
do  100  i®l , 1000 

read(18, 10,end=9999)  linein 
10  format (80A1) 

do  50  j®l,80 
lineout ( j )=linein( j ) 

if  ((linein(j) .eq."n") .and. (linein(j+l) .eq."c"))  then 
c 

c  looking  for  nccre  ncvp  or  ncendf 
c 

if  (mcont)  then 

c 

c  if  we  are  in  define  mode,  we  are  looking  only  for  ncendf 
c 

if  ((linein(j+2)  .eq.''e")  .and.  (linein(j+3)  .eq."n")  .and. 
x  (linein(j+4) .eq. "d") .and. (linein(j+5) .eq. "f "))  then 


mcont®. false. 


lineout(l)*"c" 
end  if 
else 

if  ((linein(j+2) .eq. "c") .and. (linein(j+3) .eq. "r") . and. 
x  (linein(j+4) .eq. "e"))  then 

c 

c  we  have  found  the  entrance  to  define  mode,  start  commenting  out 
c  the  lines 

c 

mcont^.true. 

endif 

if  ((linein(j+2) .eq."v") -and. (linein(j+3) .eq."p"))  then 
lineout(j+3)="g" 
endif 
endif 
endif 

50  continue 

if  (mcont)  lineout (l)="c" 
write(17,10)  (lineout(ii) ,ii=l,80) 

100  continue 
goto  9999 
9990  continue 
9999  continue 

if  (mcont)  then 
c 

c  error,  should  not  terminate  still  in  define  mode 
c 


write(6, 11) 


11 


format ( 'terminated  in  define  mode,  error  in  input  file’) 

endif 

stop 

end 


c 


8.3  Design  Constraints 

If  t  he  net  CDF  jacket  s  overhead  exceeds  20  percent  we  may  wish  to  reconsider  our  design. 
Ideally  we  hope  not  to  exceed  5  percent  on  average. 

9  Validation  Criteria 

This  section  discusses  how  we  will  know  that  the  netCDF  data  dict  ionary  and  utilit  ies  and 
subroutines  are  performing  correctly.  It  will  be  the  basis  for  the  system  tests. 

9.1  Performance  Bounds 

The  netCDF  interfaces  must  add  no  more  than  twenty  percent  overhead  to  the  model  run. 
and  we  would  like  to  see  uo  more  than  live  percent. 

An  associated  software  report  routine  will  be  called  to  document  each  netCDF  activity 
to  disk.  This  report  will  document  that  the  communications  that  we  expected  to  take  place 
did  in  fact  occur. 

9.2  Test  Cases 

Our  timing  test  case  will  run  the  models  with  ami  without  the  netCDF  interface  and  lime 
the  run  costs  of  each  to  determine  the  percentage  of  overhead.  This  will  be  done  averaged 
over  live  separate  overall  tests  to  average  out  system  load  affects,  and  over  live  separate  cases 
for  each  module  interlace. 


The  software  results  should  he  equivalent  in  all  respects,  except  performance,  up  to 
level  of  precision  used  by  the  netCDF  file  format. 


fib 
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